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Abstract 

We study magnetic vortex-like solutions lying on the spherical surface. The simplest cylin- 
drically symmetric vortex presents two cores (instead of one, like in open surfaces) with same 
charge, so repealing each other. However, the net vorticity is computed to vanish in accordance 
with Gauss theorem. We also address the problem of a flat plane in which a conical, a pseudo- 
spherical and a hemispherical segments were incorporated. In this case, if we allow the vortex 
to move without appreciable deformation in this support, then it is attracted by the conical 
apex and by the pseudosphere as well, while it is repealed by the hemisphere. This suggests 
that such surfaces could be viewed as pinning and depinning geometries for those excitations. 
Spherical harmonics coreless solutions are discussed within some details. 
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1 Introduction and Motivation 



Spheric-like shapes abound Nature: From micro to macroscale, particles and bodies with this geom- 
etry fill the Universe. Besides, spherical symmetry is very important in Physics and other sciences 
and is behind a number of triumphs of theoretical physics, namely the analytical resolution of the 
Hydrogen atom, the Schwarzschild and Kerr solutions of the Einstein field equations of General Rel- 
ativity, and so forth. 

Recent advances in Nanoscience and Nanotechnology have led, in turn, to several structures at 
extremely small scales displaying a sort of shapes: cylindrical and toroidal nanotubes, nanocones, 
and spherical nanoshells are some examples. Namely, the latter has received a great deal of efforts in 
the last years. Actually, a novel paradigm in bottom-up design for fabricating new materials lies on 
the geometry and topology of the sphere [T] , the study of optical and electronic properties of metallic 
nanoshells has revealed a number of results which could be important for spectroscopy based upon 
surface plasmon resonance [2], and spherical cavities may provide a smooth mechanism for pinning 
superconducting line vortices [3]. Concerning topological excitations in this framework, we may 
quote the formation and dynamics of vortex-like magnetisation observed in chains of spheres |1] and 
in hollowed spherical objects [5]. The relevance of topological objects lying on the spherical geome- 
try has also been considered in nematic shells textures, where a transition from half-disclinations to 
half-hedgehogs excitations is predicted as the thickness of the shells increases[6]. 

Besides these examples, excitations with non-trivial topology have a long history in Physics and 
are related to important phenomena. A classical example is the (topological) phase transition pro- 
vided by vortex-pair dissociation in a number of planar-like samples, like magnets and superfluids [7] . 
More recently, the manipulation of curling magnetisation in bi-hollowed magnetic nanodots has been 
investigated for producing basic logical operations by means of the switching of the vortex between 
the two cavities (antidots)[5]. Other pseudoparticles, generically called solitons, appear as classical 
and topologically stable solutions carrying finite size and mass in a number of non-linear models. 
For example, two-dimensional (2D) isotropic Heisenberg ferromagnets exhibit, in their continuum 
limit, spin textures of the solitonic type[9]. Such excitations are also present in theories of more 
fundamental interactions, like the skyrmions participating in the effective strong coupling [10] and 
exotic particles carrying magnetic charge[TT]. 

Solitons, vortices and other topological excitations, particularly associated with 2D magnetic 
systems, are often investigated in the lattice and continuum frameworks of Heisenberg-like models. 
Although there is a extensive literature dealing with such a subject in the planar case, the same 
cannot be said about more complex geometries, say, curved, torsed, and so on. Motivated, in part, 
by the fact that the fabrication and manipulation of magnetic surfaces with several shapes are cur- 
rently facilitated, a number of works has been devoted to such issues. For example, solitons have 
been considered in spins systems displaying the shapes of cylinders [12], [13], spheres [T3]. and cones 
|14j . Magnetic vortices of the XY model have been studied on the surface of a cone [15], in which one 
realizes that the vortex is attracted towards the apex. On the other hand, on the pseudosphere (a 
hyperbolic surface with constant Gaussian curvature) soliton appear to be unstable, decaying against 
the ground-state, while vortex-like solutions present the interesting feature of having finite energy at 
the thermodynamical limit [16]. The main lesson left by these works is that such excitations present 
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some characteristics that strongly depend on the geometrical parameters associated to the supports. 

Here, we shall be interested in analysing the basic properties of magnetic vortices whenever they 
lie on a spherical surface. It should be stressed that some analogues of this system are very important 
and have received much attention, namely in Physics of Fluids. For instance, a number of features 
accompanying the structure and dynamics of actual hydrodynamic curling fields (hurricanes, etc) in 
the Earth atmosphere are somewhat obtained from the study of vortex-like solutions on a sphere[T7]. 
[Similar plans are claimed to work for other planetary "atmospheres"]. However, it should be em- 
phasised that the latter are different from those associated to the spin field in the sense that these 
are not vectors in the physical space, so they do not experience the curvature effects directly, say, 
by means of the connection (Christoffell symbols) in the covariant derivative. Indeed, although this 
coupling is nonvanishing it does not contribute to the magnetic energy of those excitations: It is an 
extra energy called magnetoelastic[T5] which is important for studying deformable surfaces like those 
appearing in soft condensed matter problems. 

In order to carry out our investigation we shall employ the continuum limit of the Planar Rotator 
Model (0(2) isotropic model, or alternatively, the XY regime with spins confined to the internal 
plane). Among other, we realize the appearance of two singular cores for the simplest solution, so 
that the net vorticity on the whole sphere vanishes, in contrast to what is observed in the plane 
(or in other non-compact surfaces). Indeed, whenever measured isolated, the cores present the same 
charge, in such a way that they repeal each other. The energy is computed and may be written in a 
form which resembles that of a vortex in the plane, provided that the sphere radius is large enough. 
In addition, we consider a very slow (adiabatic) dynamic for a vortex on a large plane in which were 
inserted a hemisphere, a truncated cone and a part of a pseudosphere. Now, our intention is to analyse 
how these inserted geometries affect vortex motion, say, attracting or repelling it in comparison with 
the flat surface. We have seen that if the interfaces between these geometries are smooth and the 
in-plane vortex does not change its profile (no large deformation neither it develops out-of-surface 
component) then the hemisphere repeal the vortex, while in the other cases it is attracted. Attention 
is also paid to other solutions which present no singular cores, the well-known spherical harmonics, 
whose energies are analytically evaluated and are finite. Finally, we put forward our Conclusions and 
Prospects for future investigation. 

2 The continuum Heisenberg model on the sphere 

The anisotropic Heisenberg model, for nearest-neighbour interacting spins, on a two-dimensional 
lattice, is given by the Hamiltonian below: 

Hutt = -J' Yl = E + ^'^i + (1 + ^) ^tS]) , (1) 

<i,j> <i,j> 

where J' is the exchange coupling between the spins and Si = {Sf,Sf,Si) is the spin operator at 
site i. This Hamiltonian describes a number of ferro/antiferromagnets according to J < or J > 0, 
respectively. The parameter A accounts for the anisotropy: for A > spins tend to align along the 
2;-axis (easy-axis regime); for A = we have the isotropic case; while for — 1 < A < one gets the 
easy-plane regime. Finally, A = — 1 yields the so-called XY model (or the Planar Rotator Model, 
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PRM, if we focus only on 2-component spin, imposing 5*^ = 0). 



In the continuum approach of spatial and spin variables, valid at sufficiently large wavelength and 
low temperature, the time- independent model above may be written as (for that, we have considered 
a square lattice; also J = J'/2): 



= J ! ! {l + 5azX){DS''fdn (2) 



where Vt is the surface with curvilinear coordinates r/i and 772, so that dVt = \/\g\driidri2, 5ab is 
the Kronecker symbol, D is the covariant derivative, \/\g\ = ^/\det[gij]\, and gij and hat are 
the elements of the surface and of the spin space metrics, respectively, while g^'^gjk = ^l- Now, 
S = {Sx, Sy, Sz) = (sin(6) cos($); sin(6) sin($); cos(6)) is the classical spin vector field valued on a 
unity sphere (internal space), so that 6 = 6(771,772) and $ = $(?7i,?72)- The Hamiltonian above may 
be also viewed as an anisotropic non-linear a model (NLaM), lying on an arbitrary two-dimensional 
geometry, so that our considerations could have some relevance to other branches like Theoretical 
High Energy Physics and Cosmology. 

Here, our interest is to study the model above on the spherical geometry (the simplest positively 
curved surface), so that we shall be dealing with a magnetic shell approximated by a continuum 
distribution of spin sites. Thus, both the internal (spin) and the physical spaces are identical: the 
spher^. We have already parametrised the spin space by 6 and $. We also choose i!) and as the 
polar and azimuthal angles of the magnetic sphere. Then, from the relation x^ + y^ + z^ = we have 
X = R sm{'d) cos (ip), y = i?sin('^) sin(</)) and z = Rcos{'d), so that 6 = Q{'d,ip) and $ = ^{'d,ip). 
Once {R = constant, t?, (f) defines an orthogonal coordinate system, then hab and gij may be written 
as diagonal matrices; in addition, since S = {S^, Sy, S^) we get hat = Sab while the metric elements 
for the physical sphere reads: g^n) = g^^ = R^, g^^ = R^ sin^ {'d) , and g^^ = g^^ = 0. Thus, the 
Hamiltonian may be written as: 

i^sphere = J j d^ j dd [(1 + A siu^ 9) {d^Qf + siu^ 9 [d^^f] + 

^ [(1 + A sin^ 9) {d^Qf + sin^ 9 [d^^f] 

p27v pit 

J dip sm{^)d^ { [(1 + ^ sin^ ©) (^i?©)^ + sin^ © [d^^f] + 
Jo Jo 

+ [(1 + A sin^ 9) {d^ef + sin' 9 {d^^f] } , (3) 



^ Later, whenever we shall be dealing with the vortex-like solutions associated to the PRM, we shall have the internal 
spin space given by a circle, 3^ + Sy = — 1. 
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where 9^- accounts for d/drji. From this Hamiltonian there follow the (static) equations describing 
the spatial distribution and dynamics of the spin fields, like below: 



(1 + A sin^ e) 



9^(sin^9^e) + g 



—A sin 6 cos 6 
+ sin O cos O 



sini9(9^e)^ 



sin-i? 

(Ml! 

sin?? 



+ 



d^{sm^ sin^ 69^$) + d^{^ sin^ 69^$) = . 



(4) 



(5) 



As expected, the general anisotropic regime of the Heisenberg model is describe by highly non- 
linear differential equations. Suitable non-trivial solutions can be obtained provided some conditions 
are imposed. However, before proceeding further with such an analysis we should note that the 
expressions ([SllS]) resemble those counterparts for the planar (zero curvature) and for the pseudo- 
spherical (constant negative curvature) surfaces. Indeed, whenever Rsm-d is identified with r or i?r, 
while ip keeps it role as the azimuth-like angle, these expressions recover their planar or pseudo- 
spherical analogues [TB] (here, r = |?^ is the usual radial distance while Rt accounts for the distance 
measured along a pseudospherical geodesic, say, a hyperbole). Such a similarity is tied to the fact 
that all of these three geometries present constant Gaussian curvature everywhere and the measure- 
ment of their respective azimuthal angles coincides. 



To our knowledge, the isotropic Heisenberg model on the sphere was previously considered in the 
work of Ref.[T3]. There, the authors were basically interested in cylindrically symmetric solitonic 
solutions, so that A = 0, 9 = 6(^9) and $ = ^(v?) = ^ + constant. In this case the Hamiltonian ([3]) 
is simplified to: 



soliton 



2nJ 



/deV sin^e 

\1m) ^ sin^^? 



while eqs. (Ill)-(l5l) are reduced to the unique one: 

. ,d^e dQ 
smi}—— + cost}— 
dv^ dv 



sin dd - 



(6) 



sin 6 cos 
sin^9 



(7) 



whose simplest non-trivial solutions are the "hedgehog" given by G5 = ±^9 and Qas = vt ± ^9, which 
have unity topological charge {\Q\ = 1) and present finite energy equal to Sir J (we refer the reader 
to Ref . [13] for further details) . 



3 Vortex-like solutions 

Vortices are currently observed in ferromagnetic nanodisks exhibiting an in-plane structure every- 
where except at its core where the spins develop out-of-plane component (see Ref. [8] and related 
references cited therein). Whenever modelling the vortex as a continuum of spins we intend to de- 
scribe only the outer region of the core, once inside the analytical treatment is expected to give only 
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an estimative of its energy, shedding no light about its real structure and spins arrangement which 
require numeric/simulation techniques. Here, we are interested in vortex-like excitations and how 
their properties are affected by the spherical geometry. We also wish to mainly consider the so-called 
'in-plane' solutions associated to the Planar Rotator Model, whose continuum Hamiltonian may be 
obtained from ([3]) with A = — 1 and = 7r/2, and reads (once we are dealing with static solutions, 
our results apply equally well to the XY model): 



-f^PRM — J 




sm V 



sm.{'d) dd dip . 



which yields the Laplacian equation (for the region outside the cores): 

V2$(i9,^) = 9^(sin^99^$) + ^ = 0. (9) 

sm w 

If $ = ^{d^ip) then eq. above, augmented by the usual radial radial term, is solved by the well- 
known spherical harmonics which present no cores (see Section 5 for further details). The simplest 
core solution is obtained if we demand that $ has cylindrical symmetry, say, $ = $ ((/:>). In this case, 
eq. ([9]) is easily solved by: 

$((^) = g(^ + (^o, (10) 

where Q is the charge of the vortex (more precisely, of one of its cores, like below) while ipo is a 
constant accounting for its global profile, giving no contribution to its energy. The charge (vorticity) 
is formally defined in the continuum limit as: 

Q = ^ /(V$). cif, (11) 
Jc 

where the integration is evaluated along a closed path, C, around the core(s). Taking the solution 
above to the Hamiltonian ([8]) we obtain its energy: 

= 27r JQ' In ( \ = ^^^JQ'' Hcoi{§j2)) , (12) 

where = a/i? is a cutoff introduced to prevent spurious divergences at the cores, brought about 
by the continuum approach. There, a and R are the core and sphere radii, respectively. 



Note that this vortex present two cores centred at antipodal points on the sphere, say at = 
and § = 71. The appearance of the second core is not so surprising and may be understood, for 
instance, if we recall that the sphere may be stereographycally mapped to the (infinite) plane, and 
vice- versa. Indeed, on the plane the vortex energy reads (in our units) -Ev-piane = 27r J(5^(lnL — Ina), 
whose terms are the large and small cutoff, respectively. By wrapping the infinite plane to a sphere 
the vortex core centred at the origin is now taken to 'd = tt (south pole), while all those spins located 
at the infinite, L 00, are now identified with a unique one at "i? = (north pole). So, the border 
of the vortex at L — * 00 is mapped to an extra core on the sphere (a schematic view of a vortex on 
the sphere is present in Fig. [1]). At the cores the spins are expected to be directed perpendicularly 
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to the surface for minimising energy cost. Then, a more reahstic picture of the vortex is that O ^ 
as "i? ^ 0, TT while ^ 7r/2 outside the cores, hke we have obtained here (everywhere, $ = Qif + v^o? 
as above). In this approach, the cores energy may be analytically estimated by restricting the inte- 
gration in eq. ([6]) to the cores regions, say, from d = Q to da = a,/ R and from tt — da tt, in such 
a way that the spins roughly behave as a solitonic field there {da = a/R is the angular size of each 
core, like below). 




Figure 1: Top (left) and global (right) views of a vortex with Q = +1 on the sphere, showing the core at 
the north pole (where spins are expected to be perpendicular to the surface, for minimising energy). The 
arrows represent the spin field on this support, V$. 

Once da is very small, then eq. f|T2l) may be approximated to (cot(x) = + 0{x~'^)): 

^ An JQ^ \n{2R/a) , (13) 

so that the whole vortex may thought as the sum of two single-core vortices, each of charge Q, spread- 
ing throughout a hemisphere of radius R {R » a). Note also that the energy of one single-core 
vortex (vortex on one hemisphere) equals that presented by its planar counterpart in a disk of radius 
2R, so that the energy densities are the same in both surfaces. Therefore, the vortex energy on a disk 
is lower than on a hemisphere with the same radius. Then, if we insert a hemisphere into a plane 
we may expect that the former will 'repeal' the vortex from its surface (see next section for details). 
Note, however, that such an approximation must be considered with care since the divergence of the 
vortex energy is not related to the sphere radius, it rather lies possibly at the cores, as correctly 
stated by eq. ffT^ . 

By now we have pictured our vortex as a solution with two cores, centred at d = 0, tt, each of 
them carrying a vorticity Q. Then, we could be tempted of thinking on the vortex as an excitation 
with net charge 2Q. However, the situation is quite different: as a whole the vortex has a vanishing 
vorticity. Indeed, if we take two loops, each of them surrounding each of the cores, we would measure 
Qi = Q2 = Q- Nevertheless, the net charge is correctly evaluated by taking one loop containing 
the two cores. In this case, while one core is surrounded, say, in the +ip direction the another is 
circulated in the opposite sense, so that the total vorticity is -l-Q-l- {—Q) = 0. This fact is guaranteed 
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by Gauss theorem applied to the spher^. This implies that magnetic vortices lying on a spherical 
support must only appear "in pair" with the cores centred at antipodal points. It should be stressed 
that the result above does not conflit with our previous analysis of vortices and their energies, once 
we were considering the charge of each isolated core. Indeed, eq. f|T2l) gives the energy of the vortex 
whose cores are separated by s = ttR, which is the most energetically favourable configuration, like 
follows. For that, let us consider that the two cores are apart d = (for simplicity, along a circle 
with ip = constant; if they were aligned along different ip^s the case may be solved in the same lines). 
Its energy may be analytically obtained to bejf] 

(tan ( ) \ 
, , (14) 

tan(?yd/2) J 

where E^i and E^2 are the energies of the isolated cores, so that Qi = Q2, while the remaining 
term accounts for the effective potential, Vcs, between them. For preventing spurious divergences we 
take the angular separation between the isolated cores as 'dd G [da, — ^a]- Since Vefr monotonically 
decreases with increasing in the range above, then the lowest energy configuration demands the 
cores to be centred at antipodal points, say, separated by ^9 = vr (see Fig. [2]). Multi- vortex solutions 
may be considered in the same way. In this case, we can have cores with opposite charges so that 
these do attract each other, making the scenario more interesting. In the simplest case, we have two 
pairs (Q+ = +1 and Q- = —1) very close, say, around the north and south poles. An interesting 
issue to be studied is the thermodynamics involved in this dissociation of pairs and how this could 
change the spin-spin correlation function in this geometry. 



4 Slow vortex dynamics on a surface with interfaces 

For the sake of completeness, we would like to consider the problem of how a hemisphere may affects 
the dynamics of a magnetic vortex. For comparison, we also consider the effect of a cups-like defect of 
conical shape and of a pseudospherical segment (the situation is schematically shown in FiglS]). For 
that, we consider a large magnetic plane in which these geometries, with base radii R, are introduced. 
For simplifying the analysis we suppose that: i) the interfaces between the plane and each segment 
is smooth enough, so that their effects are not appreciable; ii) the inserted segments are sufficiently 
apart, exerting no influence each other, and; iii) the vortex motion is slow (adiabatic) and induces 
no changes in its cylindrically symmetric and planar profiles, like deformation near interfaces, de- 
velopment of out-of-plane component, etc. At some extent, our analysis is an attempt to carry out 
the work of Ref . [20] , where solitonic-like excitations in an line with curved sectors was considered, 
to two dimensions. 



^This result may be generalised to any vortex- like solution [T7|, and possibly for other compact surfaces topologically 
(homeomorphically) equivalent to the sphere. In the latter case since the surfaces may be geometrically distinct, further 
attention should be given to peculiar regions/points, like those where curvature blows up, for instance the vertices of 
a cube, etc. 

^Instead of explicitly evaluating an integral similar to eq. (O, we may use the fact that V<1> is analytic everywhere, 
except at the vortices cores, around which $ is a multivalued function. In this method, the result of eq. ([8]) is much 
easier and elegantly obtained. For details, see Ref. [19]. 
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Figure 2: The effective potential, V^ff/27rJ, between the two cores with charges Q = 1 on the sphere. Since 
VeS is minimised as 'd^ ~^ ^r, then the cores are located at antipodal points. 
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Figure 3: A schematic view of the segments in the plane. 



With these assumptions, only the vortex energy outside the core must be computed for our 
purposes. First, suppose that the vortex core is placed at the conical apex. In this case, its total 
energy reads (plane' in the subscript stands for the plane minus the disc area, ttE?): 

Ei= I {V^fdn+ I {V^fdxdy = 2nJQ[pln{R/pa) + lnL/R], (15) 

J cone J plane' 

which is less than -Epiane = 2nJQ'^\n{L/a), whenever /? < 1 (see Fig. Hj). Above, R/P is the size 
of the cone and (3 = sin a with a being the cone half-aperture angle. Therefore, it is energetically 
favourable to find the vortex core pinned at the apex (this attraction was previously studied in the 
work of Ref.[T5], where further details upon the conical geometry and parameters may be obtained). 
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Figure 4: The energy of the vortex (divided by 2ttJQ) as function of /3. Note that, for every /? < 1 its value 
on the cone (soUd curve) is lower than in the flat disc plane (dashed line). 



Now, let us see what is the effect of the hemisphere. In this case, the vortex energy reads: 
E2= f {V^fdn + I {V^fdxdy = 2tiJQ [ln(2i?/a) + ln(L/i?)] . 

J hemisphere J plane' 

Now, E2 > -Epiane implying that the vortex is repealed from the hemisphere. 



(16) 



Let us consider now a pseudospherical segment, of base radius R, inserted in the plane. In this 
case, the distance measured along a hyperbole from the origin to the interface with the plane is 
given hj Sr = r'TR/2 where r' is a parameter introduced for ensuring length dimension (Poincare 
disc radius) while = arcsin(2i?/r'). Therefore, the energy of the vortex if it were centred at the 
pseudospherical origin is given byo 



pseudosphere 



{V<^Ydxdy = 2nJQ 



plane' 



In 



/ tanh[arcsin(i?/r')] 
\ tanh(a/r') 



+ HL/R) 



(17) 



where a is the actual vortex core size on this geometry. Similarly to what happens in the conical 
case, E3 < -Epiane, SO that the vortex tends to be attracted and centred at the pseudospherical origin 
(see Fig. [5]). In conclusion, vortices in a magnetic plane, with hemispherical, pseudospherical and 
conical deformations sufficiently apart, are expected to be centred at the apex or at the saddle point 
of the pseudosphere (whether they actually choose the apex or the pseudosphere depends further on 



^For determining Tfj we take the projection of the pseudosphere segment on the xj/-plane. This is given by 
x"^ +y'^ = where x — r' cos{ip) sinh(T) and y ~ r' sin((p) sinh(T) while ip and r are the pseudospherical coordinates. 
In general, distances measured along the geodesies (hyperbola) are given by s = r'T/2. For further details upon this 
geometry see Ref . [T6] . 
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plane disc 

pseudosphere 




Figure 5: The energy of the vortex (divided by 2ttJQ) as function of R. Note that its value on the 
pseudosphere (solid line) is lower than in the flat disc (dashed curve), namely for R large enough. 

(5 and R parameters; such an analysis may be easily carried out). Although obtained from a very 
simple analysis it is worthy to mention that our results qualitatively agree which those observed with 
vortex lines in helium superfluids|21j. 



5 Coreless solutions: spherical harmonics spin fields 

In Section 3 we have studied the simplest cylindrically symmetric solutions of the equations of 
motion. This led us to realize that on a sphere a vortex presents two cores with same charge at 
antipodal points while the net vorticity identically vanishes. Now, if we relax the demand for such a 
symmetry and search for $ = $('(9, (/?) then we obtain the well-known spherical harmonics, say (for 
that we must start off from the complete 3D Laplacian, with radial dependence, and after all we take 
r = R = constant): 

$(i?, ^) = YJ^i^, if) = NP;;^{^) e'""^ , (18) 

where = (— 1)™-^^^^|^^^^ and P™ are the associated Legendre polynomials, defined in terms of 

the Legendre polynomials, Pnix), by P^{x) = (1 — x^)"^/^ -^Pn{x), with x = cos('i9). In addition, 
—n <m< +n and to guarantee single-valuedness of in m is integer. 

Clearly, V x {VY^) = everywhere so that no singular cores are now present. Some plots of the 
vector field VY^ are displayed in Fig. [6l Although m is important for its profile it is immaterial for 
the associated energy, like below: 

En = J J (V$)^rfv9 sm{^)d^ = 271 JN^ sm{^)d^ 



fj_pin\ I ( prn\2 
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Figure 6: Schematic views of the two simplest spherical harmonics coreless spin field solutions. On the left 
we have while on the right is showed. 



where = P™(cos('t!^)). Despite of its form this integral is easily solved for yieldingly] 

Er, = n{n^\)J. (19) 

Whether such regular structures appear and are stable on a magnetic sphere remains to be further 
investigated, for instance, by simulation and experiments. At principle, a point favouring their for- 
mation seems to be their finite energy, similar to what happens to topologically stable soliton-like 
excitations. On the other hand, energetics tends to prefer the more fundamental solutions, say, with 
small n values. 



Conclusions and Prospects 

The simplest magnetic vortex-like solution of the Planar Rotator Model, on a spherical surface, 
present two singular cores, with the same charge, which repeal each other so that they are located at 
antipodal points, for instance, at the north and south poles. In addition, its associated (finite) energy 
does not depend upon the sphere radius so that it does not blow up as the thermodynamic limit is 
taken, R ^ oo. These points should be contrasted with the common solution displaying an unique 
core, like those observed in open surfaces, whose energy diverges with the system size according ln(L). 

Furthermore, we consider a problem of slow vortex motion on a plane with inserted segments of 
conical, hemispherical and pseudospherical shapes. We realise that while the excitation is attracted 
by the conical apex and by the centre (origin) of the pseudosphere it is repealed by the hemisphere, 
suggesting these surfaces as pinning and depinning geometries for magnetic vortices. Attention was 



^We have used the identity 
/ sin(i?)di? 



2n[n + 1) (n + m)! 
271 + 1 (n-m)! 
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also paid to everywhere regular (coreless) solutions which are the well-known spherical harmonics. 
Their energies were analytically computed and shown to be finite, resembling at some extent, those 
associated to soliton-like excitations. Whether such solutions could be formed and are stable is dis- 
cussed in some details. 

As prospects for future investigation we may quote the analysis of multi-vortex solution, mainly 
in connection with a possible topological phase transition associated to the vortex-antivortex dis- 
sociation in this geometry [23j. Similar issue has been considered in the conical case, where the 
temperature transition seems to be shifted below as the cone aperture is decreased {f3 is lowered) 
[15] and on the pseudospherical support where vortex dissociation seems to take place at any arbi- 
trary small temperature [16]. The study of vortices in other geometries, like the torus, could also 
shed some light to the problem of phase transition. For instance, it has been recently verified in 
simulations that the two-dimensional Ising model defined on a toroidal lattice present two critical 
temperatures [21] instead of one, like in the usual planar case. The effects of holes (impurities/defects) 
on topological excitations lying on the spherical surface is another interesting subject to be worked 
out. For instance, hollowed nanospheres with vortex-like magnetisation could be important for po- 
tential applications, namely if they were considered at the nanoscale. 



Acknowledgements 

The authors are grateful to A.R. Moura and S.G. Alves for computational help. L.A.S. M61 and 
N.M. Oliveira-Neto are acknowledged for fruitful discussions. Professor Gary Williams (UCLA) is 
also acknowledged for having drawn our attention to some important references. They also thank 
CNPq and FAPEMIG for the financial support. 



References 

[1] D.R. Nelson, Nano Lett. 2, 1125 (2002); 
R.D. Kamien, Science 299, 1671 (2003); 

A.R. Bausch, M.J. Bowick, A. Cacciuto, A.D. Dinsmore, M.F. Hsu, D.R. Nelson, M.G. Niko- 
laides, A. Travesset, and D.A. Weitz, Science 299, 1716 (2003). 

[2] E. Prodan, C. Radloff, N.J. Halas, and P. Nordlander, Science 302, 419 (2003); 
E. Prodan, P. Nordlander, and N.J. Halas, Nano Lett. 3, 1411 (2003); 
E. Prodan and P. Nordlander, Nano Lett. 3, 543 (2003). 

[3] E.T. Filby, A.A. Zhukov, P.A.J, de Groot, M.A. Ghanem, P.N. Bartlett, and V.V. Metlushko, 
App. Phys. Lett. 89, 092503 (2006). 

[4] P. Barpanda, T. Kasama, R.E. Durnin-Borkowski, M.R. Scheinfein, A.S. Arrott, J. App. Phys. 
99, 08G103 (2006). 

[5] D. Goll, A.E. Berkowitz, and H.N. Bertram, Phys. Rev. B70, 184432 (2004). 



13 



[6] V. Vitelli and D.R. Nelson, Phys. Rev. E74, 021711 (2006). See also related references cited 
therein. 

[7] V.L. Berezinskii, Sov. Phys. JETP 32 493 (1970); ibid 34, 610 (1972); 
J.M. Kosterlitz and D.J. Thouless, J. Phys. C6, 1181 (1973). 

[8] M. Rahm, J. Stahl, W. Wegscheider, and D. Weiss, App. Phys. Lett. 85, 1553 (2004); 
M. Rahm, J. Stahl, and D. Weiss, App. Phys. Lett. 87, 182107 (2005); 

A.R. Pereira, A.R. Moura, W.A. Moura-Melo, D.F. Carneiro, S.A. Leonel, and P.Z. Coura, J. 
Appl. Phys. 101, 034310 (2007). 

[9] A. A. Belavin and A.M. Polyakov, JETP Lett. 22, 245 (1975), see also R. Shankar, J. Phys. 
(Paris) 38, 1405 (1977). 

[10] T.H.R. Skyrme, Nucl. Phys. 31, 556 (1962). 

[11] G. 't Hooft, Nucl. Phys. B79, 276 (1974); 

A.M. Polyakov, JETP Lett. 20, 194 (1974); Soviet Physics JETP 41, 988 (1976). 

[12] S. Villain-Guillot, R. Dandoloff, and A. Saxena, Phys. Lett. A188, 343 (1994); 

R. Dandoloff, S. Villain-Guillot, A. Saxena, and A.R. Bishop, Phys. Rev. Lett. 74, 813 (1995); 
R. Dandoloff and A. Saxena, Eur. Phys. J. B29, 265 (2002); 

E.A. Silva and A.R. Pereira, Phys. Status Sohdi B213, 481 (1999); Sohd State Commun. 113, 
699 (2000). 

[13] S.Villain-Guillot, R. Dandoloff, A. Saxena, and A.R. Bishop, Phys. Rev. B52, 6712 (1995). 

[14] A.R. Pereira, J. Magn. Mag. Mat. 285, 60 (2005); 

W.A. Freitas, W.A. Moura-Melo, and A.R. Pereira, Phys. Lett. A336, 412 (2005); 

V.B. Bezerra, C. Romero, and S. Chervon, Int. J. Mod. Phys. D14, 1927 (2005); 

A. Saxena and R. Dandoloff, Phys. Rev. B66, 104414 (2002), for the case of a truncated cone. 

[15] W.A Moura-Melo, A.R. Pereira, L.A.S. M61, and A.S.T. Pires, Phys. Lett. A360, 472 (2007). 

[16] L.R.A. Belo, N.M. Oliveira-Neto, W.A. Moura-Melo, A.R. Pereira, and E. Ercolessi, ^^Heisenberg 
model on a space with negative curvature: Topological spin textures on the pseudosphere" , Phys. 
Lett. A (2007) [in press], doi:10.1016/j.physleta.2007.01.044 . 

[17] D. Crowdy and M. Cloke, Phys. Fluids 15, 22 (2003); 
D. Crowdy, J. Fluid Mech. 498, 381 (2004). 

[18] A. Saxena, R. Dandoloff, and T. Lookman, Physica A261, 13 (1998). 

[19] P.M. Chaikin and T.C. Lubenski ^^Principles of Condensed Matter Physics" , Cambridge Univ. 
Press; 1st edition (1995), see mainly Section 9.3; 

G.A. Willians and E. Varoquaux, J. Low Temp. Phys. 113, 405 (1998). 
[20] R. Dandoloff and A. Saxena, Phys. Lett. A358, 421 (2006). 



14 



[21] P. Voll, N. apRoberts- Warren, and R.J. Zieve, ''Surface curvature and vortex stabilit'i/', 
|cond-mat/0601395| 

V. Vitelli and A.M. Turner, Phys. Rev. Lett. 93, 215301 (2004). 

[22] G.B. Arfken and H.J. Weber, ''Mathematical Methods for Physicists" , sixth edition, Elsevier 
Academic Press, 2005; see Chapter 12, mainly pages 771-802. 

[23] V. Kotsubo and G.A. Wilhams, Phys. Rev. B33, 6106 (1986); 
H. Cho and G.A. Wilhams, J. Low Temp. Phys. 110, 533 (1998); 
G.A. Williams, Phys. Rev. B73, 214531 (2006); 
W.A. Moura-Melo, work in progress. 

[24] I. Hasegawa, Y. Sakaniwaa and H. Shima, J. Mag. Magn. Mat. 310 (2007) 1407. 



15 



